ON THE RAYLEIGH-TAYLOR INSTABILITY FOR THE 
TWO-PHASE NAVIER-STOKES EQUATIONS 
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Abstract. The two-phase free boundary problem with surface tension and 
downforce gravity for the Navier-Stokes system is considered in a situation 
where the initial interface is close to equilibrium. The boundary symbol of 
this problem admits zeros in the unstable halfplane in case the heavy fluid 
is on top of the light one, which leads to the well-known Rayleigh- Taylor 
instability. Instability is proved rigorously in an Lp-setting by means of an 
abstract instability result due to Henry |12| . 



1. Introduction 

Of concern is the motion of two immiscible, viscous, incompressible capillary 
fluids, fluidi and fluid2, that occupy the regions 

n,{t) = {(x, ?/) e R" X R : {-iy{y - h{t, x)) > 0, i > 0}, i ^ 1, 2. 

The fluids are separated by a sharp interface 

T{t) = {{x,y) xR: y ^ h{t, x), t > 0} 

with an unknown function h that needs to be determined as part of the problem. 
The motion of the fluids is governed by the incompressible Navier-Stokes equations 
with surface tension and downforce gravity and reads as follows, where i — 1,2] 
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The constants pi > and pi > denote the densities and the viscosities of the 
respective fluids, a stands for the surface tension and 7a is the acceleration of 
gravity. Moreover, S{u^ q) is the stress tensor defined by 

S{u,q) = fit{Vu+ (Vu)'^) - ql in n,{t), 

and |u] = — ^|s7j(t)) lr(t) denotes the jump of the quantity v, defined on the 

respective domains fli{t), across the interface T{t). Finally, k — K.{t, •) is the mean 
curvature of the free boundary T(t), v — v(t^ •) is the unit normal field on r(t), and 
V — V(t^ •) is the normal velocity of r(i). Here we use the convention that v(t^ •) 
points from into (^)i f^nd that k(x,€) is negative when r2i(t) is convex 

in a neighborhood of a: G r(t). System (jl.ip comprises the two-phase Navier- 
Stokes equations with surface tension subject to gravity. In order to economize our 
notation, we set 

P = PiXfiiit) + P2Xo2(t): M = A^iXni(t) + A*2Xo2(t)) 

where x denotes the indicator function. It is convenient to introduce the modified 
pressure q := q + p^aU- With this convention system (jl.ip can be recast as 

' p{dtu+{u-V)u) - fiAu + \7q = in n{t) 

div u = in Q{t) 

-lS{u, q)^] = anv + |p]7a?/ on T{t) 

|u]=0 on T{t) (1.2) 

V ^u-v on T{t) 

u(0) = uq in ^^0 

r(o) = To . 

Given are the initial velocity uq : ^ ]R"+^ with fig := f^i(O) U 512(0) as 
well as the initial position Fq = graph (/iq). The unknowns are the velocity field 
u(t, •) : U,{t) E"+^, the pressure field q{t, •) : n(t) R, and the free boundary 
F(t), where n{t) := Qiit) U 5^2 (i). 

In case that ^i{t) is a bounded domain, ja = 0, and 02(i) = 0, one obtains the 
one-phase Navier-Stokes equations with surface tension, describing the motion of 
an isolated volume of fluid. For an overview of the existing literature in this case 
we refer to the recent publications [17l [HI [20l [21] . 

The motion of a layer of viscous, incompressible fluid in an ocean of infinite extent, 
bounded below by a solid surface and above by a free surface which includes the 
effects of surface tension and gravity (in which case flo is a strip, bounded above 
by Fq and below by a fixed surface F;,) has been considered by jTll^ [51 [M] . 

If the initial state and the initial velocity are close to equilibrium, global existence 
of solutions is proved in 2 for cr > 0, and in 24J for cr > 0, and the asymptotic 
decay rate for i — > cx3 is studied in We also refer to ^Sj, where in addition 
the presence of a surfactant on the free boundary and in one of the bulk phases is 
considered but gravity is neglected. 
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Previous results concerning the two-phase problem (II. 2|) with 7q = in the SD- 
case are obtained in [71 13 [9l |22] . In more detail, Densiova [8] establishes existence 
and uniqueness of solutions (of the transformed problem in Lagrangian coordi- 
nates) with V G 14^2 *'''^ for s G (5/2, 3) in case that one of the domains is bounded. 
Tanaka [22] considers the two-phase Navier-Stokes equations with thermo-capillary 
convection in bounded domains, and he obtains existence and uniqueness of solu- 

s s /2 

tions with {v,9) G Wj' for s G (7/2,4), with 9 denoting the temperature. 
Here we are interested in the situation where Tq is close to a plane, say R" with 
n > 2, i.e. Tq is a graph over M", given by a function Hq that is small in an appro- 
priate norm. Then it is natural to transform the problem to a flat fixed interface, 
and solve the resulting quasilinear evolution problem. Our basic well-posedness 
and regularity result for problem (|1.2p reads as follows. 

Theorem 1.1. Fix p > n + 3 and let 

iuo,ho) G X W^3-2/p(^n) 

be given. Assume that the compatibility conditions 

divuo = on fta, l^DaVQ - ^{i^q ■ DQi^Q)voj = 0, [uo] = on Tq, (1.3) 

are satisfied, where Dq — (Vug -|- (Vito)^), and Vq is the unit normal field on Fq. 
Then for each a > there exists rj > such that for 

there exists a unique classical solution (u,q,T) of problem (jl.2p on (0, a). In 
addition, Ai = ^te{o,a){{'t} x r(f)) is a real analytic manifold, and the function 
{u,q) : O M"+^ is real analytic, where O := Uig(o.a)({t} x fl{t)) . 

Proof. This result is proved in [16j in case that 7a = 0. The proof given there 
extends to the case 7a > 0, as the additonal term |yo]7a/i on the interface is of 
lower order. Actually, in Section 6 we shall give a different existence proof based 
on the implicit function theorem. □ 

We mention that system (jl.2|) has also been analyzed in HI] for initial data that 
are not necessarily close to equilibrium. More precisely, it is proved in [17^ that 
(II. 2p admits unique solutions (on a possibly small time interval) that have the 
same regularity properties as above, provided ||V/io||oo is small enough. 

It is the purpose of this paper to prove mathematically rigorously that the trivial 
solution (u, h) = (0, 0) of problem (|1.2p is unstable in the phase manifold VM, to 
be defined below, in an ip-setting in case that the heavy fluid overlies the lighter 
one, i.e. if p2 > pi- This is the Ray leigh- Taylor instability which is well-known 
in Physics and Hydrodynamics, cf. [4, 6, [Til [IS] the references given there. 
The Rayleigh- Taylor instability manifests itself in the way that any disturbance 
of the equilibrium solution (w, h) = (0, 0) will grow to produce spikes of the heavy 
fluid moving downward and bubbles of the light fluid moving upward. The precise 
statement of our main result is as follows. 
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Theorem 1.2. Let p > n + 3. Suppose that pi, p2, /^i, /^2, cr, 7a > are constants 
and p2 > pi- Then the trivial equilibrium (u,q,h) — (0,0,0) is Lp-unstable. More 
precisely, there is a constant Eq > such that for each S > there are initial values 

subject to the compatibility conditions (|1.3p in Theorem 1.1 with 



iFol 



w; 



such that the solution (m, h) for some Iq € (0, a] satisfies 
\W{to)\\^2-2/p + \\h{to)\\^3-2/p > eo- 

Our method depends on the proof of Theorem 1.1 presented in [T6|, as well as on an 
abstract instability result for iterates of a mapping due to Henry |l2j , applied here 
to the Poincare map or time-one-map of the system. To verify the assumptions 
in Henry's result we show that the boundary symbol s(A, r) admits zeros (Aq, tq) 
in the unstable half-plane in case p2 > pi and prove that such a zero induces 
the spectral values Aq for the linearized operator of the problem at the trivial 
equilibrium. 

2. Reduction to a Flat Interface and Linearization 



The nonlinear problem (jl.2p can be transformed to a problem on a fixed domain 
by means of the transformations 

v{t, X, y) := {ui, . . . , Un){t, x,y + h{t, x)), 

w{t, X, y) := w„+i(t, x,y + h{t, x)), 

TT{t, X, y) := q{t, x,y + h{t, x)), 

where t E J = [0,a], x E M", y G K, y 7^ 0. With a slight abuse of notation 
we will in the sequel denote the transformed velocity again by u, that is, we set 
u= (v,w). With this notation we obtain the transformed problem 



pdtu — pAu + Vtt 


^ F{u,n,h) 


in 
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^Fd{u,h) 


in 




-Ipdyvj - IpV^wj 
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on 
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(2.1) 



for t > 0, where R"+^ := {{x,y) e R" x R : y ^ 0}. More details on this 
transformation, on the nonlinear right hand sides, can be found in [161 . Here we 
should point out, however, that the definition of Gy in this paper differs from that 
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in [T^ in the following way: solving the second line of formula (2.7) in [16] for 
[tt] and substituting the result into the expression for Gv in formula (2.8) of [16] 
results in 



+ {lfldyWj - iVh\ l^iV,w]) + \Vh\^^idyW]}Vh. 



(2.2) 



Thus the quantity [tt] can be eliminated in the nonlinearity G„. 

The linearization of (|2.ip at (m, h) — (0, 0) leads to the linear inhomogeneous 
problem 



pdtu — ^iAu + Vtt 
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(2.3) 



pTl + l 



(2.4) 



We are interested in the regularity class 

u e ffi(J; Lp(M"+i,M"+i)) n Lp(J;il2(M"+i, 

7reLp(J;iji(]R"+i)), 

where J — [0, a]. In the following, denote as usual the Sobolev spaces if 
TO G Z. For non-integer s, are the Sobolev-Slobodeckii spaces, and the 

Bessel-potential spaces. For K G {H,W}, by Kp we mean the homogeneous 
version of K^. Note that — Wp for integer values of s, but that in general 
these spaces are different. We refer to flFj Section 2] for more details. 

If we assume a solution in the class (|2.4p . then for the right hand sides / and 
fd we necessarily have / e Lp{J x M"+^,M"+^) and 

/, G H'piJ;Hp\W^+')) n Lp(J;Hi(i"+i)), 

since the operator div maps Lp into H^^^. By trace theory we necessarily have 



), and the lateral trace of u belongs to 



and that of djU to 



R"+i)), 



Yi Wy^-^/^P{J-Lp{R",R''+'^))nLp{J;W^-^/PiR'',M 
see for instance [10] . Therefore gv & Yi, and if in addition 

H e Lp(R")) n Lp{J; w^-'/piR")), 
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then we also have that g^, G Yi. 

Concerning the regularity of the height function h we note that the equation 
for h lives in the trace space Yq, hence naturally h should belong to 
h e VF"2-i/2p( J. n Fi( J; V72-i/p(M")). 

On the other hand, the equation for the normal component of the normal stress 
lives in Yi, and contains the term Ah, hence h should also belong to the space 
Lp{J;Wp~^^P (W^)). These considerations lead to the following natural space for 
the height function h 

This then implies gh (zYq, as well as ho € Wp^^^P{W^) by trace theory. Our next 
theorem states that in this setting, problem (j2.3p admits maximal regularity; the 
described regularities of the data are also sufficient. In particular, the solution 
map defines an isomorphism between this space of data and the solution space 
defined above. 

Theorem 2.1. Let 1 < p < oo be fixed, p ^ 3/2, 3, and assume that a, ja, Pi o,nd 
fj,i are positive constants for i — 1,2, and set J — [0,a\. Then the instationary 
Stokes problem with free boundary p.3p admits a unique solution (u,7r, |7r],/i) with 
regularity 

ue H^{J;Lp{«"+\W'+^))nLp{J;H^{R''+\W'+^)), 
7reLp(J;iji(R"+i)), 

H e iy;/2-i/2P( J; Lp(R")) n LpiJ; M^pi-i/f (M")), 

h € VFp2-i/2p( J; Lp{W')) n Hl{J; W^p2-i/p(R")) n Lp{ J; W^~^/p{K")) 

if and only if the data (/, fd, g, gn, uq, ho) satisfy the following regularity and com- 
patibility conditions: 

(a) /eLp(J;Lp(R"+i,M"+i)), 

(b) fd e H^piJ; ij-i(K"+i)) n Lp{J; ffi(R"+i)), 

(c) g = {g^ng^o) G T</'-'/'''(J;ip(R",R"+i))nLp(j;M^;-^/^(R", 

(d) gh e T<"'/'^(J; Lp(R")) n LpiJ; W^-'^^iW^)), 

(e) no e iyp"^/''(R"+\R"+i), ho G Wp^^^^iW), 

(f) divuo = /d(0) m R"+i and |uo] on R" i/p > 3/2, 

(g) -Ifidyvoj - I/iV.u'o] = 5^0) on R" i/p > 3. 

The solution map [[f,fd,g,gh,uo,ho) i~-> (u, tt, [tt], /i)] is continuous between the 
corresponding spaces. 

Proof. For a detailed proof of Theorem 2.1 in case 7q = we refer to [16]. The 
proof carries over to the case 7a > since the term |p]7a^ is of lower order. The 
only change occurs in the boundary symbol s(A, |^|); see Section 3. □ 
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We also need a corresponding result for the stationary linear problem 



/oA*w — /iAu + Vtt 
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in 
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on 
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(2.5) 



where A* > is sufficiently large. It reads as follows. 

Theorem 2.2. Let 1 < p < oo fee fixed, and assume that a, ^a, Pi o,nd fii are 
positive constants for i = 1, 2, and that A* > is large enough. Then the stationary 
Stokes problem with free boundary (|2.5p admits a unique solution (u, tt, h) with 
regularity 

if and only if the data {f , fd, 9, gh) satisfy the following regularity conditions: 

(a) fd e iyp'-'/''(E"+i) n ijp-i(E"+i), 

(b) <? = (.g.,.g^) e iyp'-'/''(R";E"+i), 

(c) <?„eT4^r'/^(R")- 

The solution map [{fd, 9, 9h) ^ [u, tt, [tt], h)] is continuous between the correspond- 
ing spaces. 

Proof. The proof will be given at the end of the next section. □ 



3. Zeros of the Boundary Symbol 
As shown in our paper [17j , the boundary symbol of the linear problem is given 

by 

s{X, \Q) . A + ^ ^^^^^1^1 ^ ^^^^^^^^^ ^ . (d.l) 

Here A denotes the co- variable of time t and ^ that of the tangential space variable 
X S M", and we employed the abreviations 

The boundary symbol s(A, |^|) has been studied in detail in the papers [1^1111] in 
case 7a = 0, and in [T^ for 7^ > 0. It has been shown in [T71 Remarks 3. 2(b), (c)] 
that s(A,t) does not admit zeros (A,t) 7^ (0,0) with ReA > and r > in case 
P2 l£ Pi, i.e. in the stable case. On the other hand we have the following result. 
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Proposition 3.1. Suppose p2 > pi > and pi, fj,2, ly^Ja > are constants. Then 
for each ^ G M" with 

< lei < n := (74pl/a)i/2 

there is A(|^|) > such that s(A(|^|), |^|) = 0. Every zero of s{\, |C|) with ReA > 
is real. 

Proof. Note first that s(A, 0) 7^ unless A = 0, hence we may assume ^ 7^ below. 
It is convenient to the use the scaling C, — X/t"^ where r = |^| G (0, 00). By a slight 
abuse of notation we set 



and obtain 

s(A,t) =r2(C + V(r)fc(C)), 



where 



and 



{Pi+P2)t {P1+P2W' 

4 m (0^2(0 



C + 



Thus (A, r) is a zero of s if and only if (C, r) satisfies C + 4'{T)k{C) = 0. It has been 
shown in that argfc(C) G (-7r/2,0] if argC S [0,7r/2]. This implies that for 
ReC > 0, C 7^ we have C,+'^{T)k{C,) ^ if either r > or C is non-real. Thus we 
need to show that for r G (0, r*) there is a a zero (^{t) > of C + ?A(T)fc(C) = 0. But 
the function ^(C) := C/^(C) is nonnegative and strictly increasing for C ^ 0, it is 
zero at C = and behaves like as t — + cx3, hence $(0, 00) D (0, 00). On the other 
hand, — ■(/'(O, t*) = (0, 00) imphes that for each r G (0, t*) there is a unique ^(r) > 
such that <I>(C(r)) = — ■(/;(t), i.e. (C(T),r) is a zero of the function C, + V'(''')fc(C)j 
which yields the zeros (C(|C))I^Pi of tt^^ boundary symbol s(A, |^|). □ 



The critical value t» = -v/IplTa/c is known as the cutoff wave number in the 
literature, cp. e.g. [IT]. Note that $(C) as C — > 00, while $(C) ^ (/k{0) as 

C ^ 0. This gives the asymptotics 



C(t) ~ Vlpha/{pi+P2)r-'^', r^O; CM ^ c(r, - r), r ^ t„ 
with c = 2(Tfc(0)/(pi + P2)t^, which implies for A(t) — t^z(t) 

A(r) ~ V[pHJW+P2)r'/^ r 0; A(r) ^(r, - r), r ^ r,. 

Thus A(t) ^ as t ^ 0, n, hence the function [r ^ A(t)] has a positive absolute 
maximum in the interval (0,r,) which we denote by Aqo > in the sequel. 

Proof of Theorem 2.2: Necessity is obtained by trace theory. To prove sufficiency, 
we can use the same reductions as in Sections 3-5, with the notable difference 
that here we only need to consider the stationary case with A* a fixed parameter. 
As in the proof of [16, Theorem 5.1] it then remains to consider the boundary 
symbol s(A*, |^|). We have seen above that for A* > Ac© the boundary symbol 
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does not vanish, and we may estimate s(A,, t) from above and below by Cj{l + |t|) 
with appropriate positive constants Cj, j = 0, 1. This is vaUd for 

reS^-lCeC: Ct^O, largCl <ry} 

with sufficiently small 77 > 0. Therefore by Mikhlin's Fourier multiplier theorem 
(1 + Id) defines a bounded linear operator in i/^(R") and in W^{W) for 

all s e M. □ 

4. The Spectrum of the Linearization 

As a base space for the functional analytic setting we use 

Xo - Lp^<,(M"+^M"+i) X Wp2-i/P(M"), 

where the subscript a means solenoidal, and we set 

Xi = i/2(R"+i;M"+i) X W^-^^P{R") 

As before we use the decomposition u = {v, w). Define a closed linear operator in 
Xq by means of 

A{u,h) = {-^Au+-Vn,-w), (4.1) 
P P 

with domain Xi := D(A) C Xi 

D{A) = {(m, /i) e Xi n ATo : M = 0, {^idyV + iiW^w\ on K"}. 



The pressure field tt in (|4.ip is determined as the solution of the transmission 
problem 

(iv7r|V0)i, = (^Azi|V0)i„ e W^,MK"+'), ^^2) 

1^1 - 2|M9yU;l + rxA/i + [ph./i on M" . 

One should compare this operator A with 20J, where the corresponding operator 
in the one-phase case has been introduced and analyzed. 
Concerning the transmission problem, we set 

Ti- := tt/p, /:=pAu/p, 5 := 2|/^9yw] + crA/i + |p]7aft.. 

Then the solution of the transmission problem consists of two parts tt = tti +7r2 . tti 
is given by tti := — div_D~^/, i.e. Vtti = Rf , where R denotes the Riesz-transform 
with symbol ^ ® C/ICP- On the other hand, 112 solves 

A7r2 = in R"+i, [^^2! =9- bvn] 50 and [9,7r2] = on R". 
The solution of the latter problem is given by 

■^2{-,y) = 7 \ 7 sign(y)F(|y|).9o, 

[Pi + P2) 

where {P{s) : s > 0} denotes the Poisson semigroup on M". Thus Vtt e Lp(M"+^) 
since q and 50 belong to Wp ^'^{W^) and / e Lp(K"+i), see for instance formula 
(2.14) in [in] for the assertion for ^2- 



10 



J. PRUSS AND G. SIMONETT 



System p.3p can be rewritten as the abstract evolution equation 

j^iu,h) + A{u,h) = {fd,gh), t>0, (w(0),/i(0)) = K,/io), (4.3) 

provided {fd, g) = (0, 0). Since (|2.3p has maximal Lp-regularity, the abstract prob- 
lem (j4.3|) has maximal Lp-regularity as well. In particular, —A generates an an- 
alytic Co-semigroup in Xq. Concerning the spectrum <t{A) of A we have the 
following result. 

Proposition 4.1. Suppose p2 > pi > and p-i, lJ-2,'^,la > are constants, and 
let A with domain D{A) — Xi be defined in Xq as above. 

Then Aq G (^{-A) for each zero (Ao,|Co|) e (0,oo) x (0,cx)) o/ s(A, |^|) = 0. In 
particular, a{~A) n 0+ = [0, Aoo] with Aoo > from the previous section. 

Proof, (i) The idea is to show that whenever we have a zero (Ao,|^o|) of the 
boundary symbol s then Aq is an approximate eigenvalue of —A, hence in a[—A). 
So suppose s(Ao, |fo|) = 0. We define functions by means of 

h,{x) = e*««-^xe(a;), X e M", e > 0, 

where x is a Schwartz-function such that its Fourier-transform x is a test function 
with supp (x) C i?Rn(0, 1), x(0) — 1 and Xe(x) '■= xi^^)- This implies by means 
of the formula h^{^) — e^"x((C ^ Co)/e) that the support of is contained in the 
ball Bk^{^o,£)- In particular, for all r G R, 1 < p < oo and K £ {H,W}, the 
operators £>„ :— —A, dI/^ and Dn^^^ as well as Bq :— XqD~^ act boundedly in 

K;^M") ■■= {g e A-(R") : supp (5) C SK„(eo,£)}, 

1 /2 

as long as, say, e < |Co|/2. Note that the spectra of Dn and Bo in ^ are 
contained in [a, /3] for some < a < (3 < oo. 

(ii) The Lp-norm of is l/ielp = £ "^^llxllp? ^^'^ moreover we have 

\\he\\H-^ < ce-"/P\\x\\H^^, k<A, 0<e< |Co|/2. 

Hence by the interpolation inequality one sees that 

coe-"/^ < < Cle-"/^ 0<e<|eo|/2, re [0,4], 

where K G {H, W} and cq, ci are positive constants. 
Next we solve the Stokes problem 

pXou - nAu -I- Vtt = in M"+i 



div u — 



m 



pri+l 



-lfidyvj-l^i\/,wj^O on R" (4.4) 

-2[/i9,u;] + H = (aA + |p]7„)/ie on R" 
H = on R". 

This is possible since Aq > 0, and we obtain for the solution (u^, tt^) the estimate 

||ue||//2(R„+i) + ||V7re||L^(R-+i) < C\\h,\\^3-i/p,^„. < Ccie-"/P. 
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We then have 

Ao(Ue, hs) + A{Ue, he) = (0, Qe), 

where = XoK + {gD\1'^ - [pha-Dn ^^^)(pi + P':^~^'k{\QD-^)he. 
(iii) It remains to estimate g^- First we observe that 

hence 

and therefore 

(C - Bo)/ie = (C - Ao|Cor')/i. - £|^or'Sor„ 

which yields 

(C - B^Y^he - (C - Mi^VY^h, + £|eor'(C - Ao|eor')-'(C - b^y^b^t,. 

With = (cr — CH7a/Ao)(/Ci + p2)~^k{Q, by means of Dunford's functional 

calculus in /^^^(R") this yields, with a closed contour F in the open right half-plane 
surrounding [a, 13], 

h{Bo)he = 2^ ^ fci(C)(C - BoY'he dC 

= hiXolioY'Ye + l^hiOiC - Xo\ioYT'iC - BoY'Bor.dC 

= ki{\o\^oY^)he+eKire, 

where the operator Ki is bounded and does not depend on e. Here we used the 
fact that fci(C) is bounded. Finally, we have in a similar way 

Dl/'K = ^ / V^(C - DnY'he dC 
= \^o\he+eK2re, 

with a bounded operator K2 that is independent of e. Thus in summary we have 
g, = (Ao + {aDi/^ - lph,D-y^){p, + p^Y'KBo))^ 
= s(Ao, + e(fci(Ao/|^o|')i^2 + KiDl/^)re 

since by assumption .s(Ao, [^o|) = 0. The operators Kj are bounded and Dn is so 
on functions whose Fourier transform has compact support bounded away from 0, 
hence we obtain 

ScaHng the functions by the factor e"/^' we have /ie)||js:o > c > and 
(Ao + A){ue, he) ^ as £ — > 0. Thus Aq must belong to the spectrum of {—A). 
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(iv) Finally, if A G C+ is not in [0, Aoo], then s(A, •) does not vanish. Therefore 
the boundary symbol can be inverted by means of Mikhlin's multiplier theorem 
and then solving the remaining Stokes problem we see that A + A is boundedly 
invertible. Thus such A belong to the resolvent set of ^. The proof is complete. □ 

5. Henry's Instability Theorem 

For the reader's convenience we provide the statement of Henry's instability the- 
orem, see jl2l Theorem 5.1.5], for a more specialized situation which is appropriate 
for the problem under consideration in this paper. 

Theorem 5.1. Let X be a Banach space, Te C^{Bx{0, g); X) such that T{0) = 
and the spectral radius spr T'(0) o/T'(0) is greater than one. Then the origin is 
unstable in the sense that there is a constant Eq > such that for each S > there 
is xs G Bx{0,S) and N E N such that the sequence Xk ■= T^xs G Bx{QtQ) is 
well-defined for a/Z fc G {1, . . . , and xs\ > eo- 

As an illustration consider as in [14 and [18 the quasi-linear evolution equation 

u + A{u)u^ F{u), t> u(0) = uo, (5.1) 

where Xi ^ Xq are densely embedded Banach spaces, X^ := {Xq, Xi)i_i/p p is a 
real interpolation space of order 1 — 1/p and power p G (l,oo) between andXi. 
The nonlinearities {F^A) : X^ — > Xq x C{Xi,Xq) are of class . Assume that 
M* G Xi is an equilibrium of (j5.ip . i.e. A(u^)u^ = F(u^,). Assume further that the 
linearization of ()5.ip has maximal Lp-regularity, i.e. the operator Aq defined by 

Aqv A{u^)v + [A'{u^)v]u^ - F'{u^)v, v £ Xi, 

is 7?.-sectorial. Due to the results in [T3] and [TB], the Poincare map of (|5.ip . 
[uq Tuq — :— u(a, uq) — u*], where u{t, uq) denotes the solution of (|5.ip with 
uo G Bx (u* , p) , is well-defined and satisfies the assumptions of Henry's instability 
theorem in X := X^, provided p > is sufficiently small. The derivative r'(0) 
equals e~^°°', which by the spectral mapping theorem for generators of analytic 
Co-semigroups has spectral radius greater than one if and only if the spectrum of 
the operator —Aq contains points in the open right half-plane. In particular, in 
contrary to [14] no spectral gap is required. Therefore we may conclude that the 
equilibrium is unstable in the natural phase space X^ of the problem. 

In the situation of the two-phase Navier-Stokes equations with surface tension and 
subject to gravity, things are more involved due to the inherent nonlinearity of the 
compatibility conditions for the transformed problem 

div u ~ Fd(u, h), 

-It^dyv] - IfiV^w] ^ Gyiu, h), (5.2) 

-2[t^dywl + H - [aAh + lplja)h = G^(u, h), 

on the interface. This leads to a nonlinear phase manifold where the semi-flow lives 
on. The way out is to parameterize this manifold. Nevertheless, the linearization 
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of the time-one map will turn out to be the operator e"^ with A described in 
the previous section, which by Proposition 4.1 has spectral values in the open 
right half-plane. This way we will still be able to apply Henry's result, and as a 
consequence to obtain a rigorous proof for the Rayleigh- Taylor instability. 



6. Proof of the main result 

(i) In a first step we parameterize the phase manifold VAi of system (|2.ip which 
lies in 

:= {{u,h) e W"2-2/p(]jn+i.Rn+i) ^ W^-^/PiR") : {u] = 0} 
and is defined by 

VM := {{u, h)<EX^: divu = w • V/i in R"+\ -l^l^yvl - = G.^{u, h)]. 

We want to parameterize VM locally near (0, 0) over the closed linear subspace 
C X^, defined by 

X° {{u, h)eX^: divM in M"+\ {fidyvj + l^lV^wl = on M"}, 

by means of an analytic map $ : Bxo (0, r) VM which is bijective onto its range 
such that $'(0) = /. For this purpose we consider the problem 

/9A*M — /iAw + Vtt = 

div u — Fd{u + u, h + h) 

— {^dyvj — l^VxwJ ~ Gi,{u + u, h + h) 

-2lfidyiw + w)j + H = (aA + IpjjaKh + h) + G,„(ii + n,h + h) 

M = o 

A,/i — w\y^Q = H{u + u,h + h) 

with some fixed, sufficiently large A, > and given (u, h) E Bxo{0, r) C We 
write this equation in short hand notation as Lx^z = N{z + z) in 

{(w,7r, H,/i) : {u,h) e X^, tt G W^^-^/f (M"+i), H G W^-'/P{Rn} 

where z ~ (u,7r, |vr],/i) and z = ({t,0, 0,/i). Since the nonlinear terms are poly- 
nomial, it is not difficult to verify that N is real analytic and N'{0) = 0; see p!6l 
Proposition 6.2] for a related result. is invertible by Theorem 2.2, and hence 
the implicit function theorem yields a unique solution z = '^{u,h) G Z-^ near 0. 
The mapping 

with r chosen small enough, is real analytic and satisfies 3'(0) — 0. Denoting by 
P : Z-y — > Xry the projection given by Pz :— (w, h) for z = (u, tt, [tt], /i), we set 

<^{u,h) := {u,h) + (t){u,h) with (i>{u,h) -.^ P'i{u,h). (6.1) 
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Then $(5^0 (0,r)) C VM, $ : Bxo{0,r) is real analytic, $'(0) = /, and $ 

is injective. Hence it remains to show local surjectivity near 0. So suppose that 
("0, h) G VM has sufficiently small norm. Solving the problem 



pri+l 



pA*u — /iAu + Vtt = in 

div u = /i) in 

-iM^y^'l - l^J■^xwl = G,(w, /i) on M" 

-2|Ai9j,i(;] + H-((7A + Ipl7„)/i = G^(jI,/i) on M" 

|u] = on M" 

A*/i - w|,j=o ^^(u, h) on M", 

by means of Theorem 2.2 yields a unique solution z — (u, tt, |7r],/i) G Z^. One 
readily verifies that {u, h) := {u — u,h — h) belongs to X° and (f){u, h) — {u, h), 
showing surjectivity of <i> near 0. In particular, VM C is a real analytic 
manifold near G X-^. We observe that due to 

-2lfidyiw + w)] + H = (ctA + lpha){h + h) + G^{u + u,h + h) on M" 
the last condition in (|5.2p is satisfied as well. 

(ii) Now we proceed as in the proof of [TBI Theorem 6.3], employing the notation 
of that proof. For a given [uq, Hq) G VA4 we construct the extension z^{uq, Hq) G 
El (a). The map 

[iuo,ho)) ^ z4$(uo,M)] : 5jfo(0,r) ^E(a) 
is real analytic. Therefore, fixing a > the mapping 
if:oE(a)xBxo(0,r)^oF(a), 

(6.2) 

K{z, {uo, ho)) := L{z + z^^uq, ho))) - N{z + z,{<^{uo, K))) 

is real analytic. We have -ftr(0, 0) = as well as DzK{Q, 0) — L. Therefore the 
implicit function theorem yields a real analytic map 

*:Bxo(0,r)^oIE(a) 

such that K{'i'{uo,hQ),{uo,ho)) — for all {uo,ho) G Bxo{0,r), with r chosen 
sufficiently small. Thus 

z := ■i'{uQ,ho) + z^{uo,ho) eE{a) (6.3) 
is the unique solution of problem (2.1) and, moreover, the mapping 

[{uo, ho) ^ z{uo, ho)] ■■ Bxo (0, r) E(a)) (6.4) 

is real real-analytic. 

(iii) Having obtained a unique solution z = z{uo,ho) G E(a) of ()2.ip . we can 
employ the same arguments as in steps (vi)-(vii) in the proof of [16, Theorem 6.3] 
to establish analyticity of the solution as stated in Theorem ll.il 
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(iv) Differentiating the mapping [{uo,ho) i"> z(uo,ho)] w.r.t. the initial value 
{uQ,ho) one sees that the linearization at (0,0) of (|2.ip is given by the Cauchy 
problem 

j^iu, h) + A{u, h) = 0, (m(0), h{fi)) - (mo, ho). 

This implies that the linearization at (0, 0) of the time-one-map [{uq, Hq)] i— > z{l)\ 
of the nonlinear problem (j2.ip is given by e^^ . Since {—A) generates an analytic 
Co-semigroup in Xq the spectral mapping theorem yields spr (e^"*) = e^°° > 1 
by Proposition 4.1, hence we may apply Theorem 5.1 to obtain instability of the 
trivial solution. The proof of our main result is therefore complete. □ 
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